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Abstract 

In a scenario in which fourth generation fermions exist, we study effects of new physics on the differential decay width, 
forward-backward asymmetry and integrated branching ratio for B — » X a £ + £~ decay with (£ — e,/i). Prediction of the new 
physics on the mentioned quantities essentially differs from the Standard Model results, in certain regions of the parameter 
space, enhancement of new physics on the above mentioned physical quantities can yield values as large as two times of the SM 
predictions, whence present limits of experimental measurements of branching ratio is spanned, contraints of the new physics 
can be extracted. For the fourth generation CKM factor V^ b V t i s we use ±10" 2 and ±10" 3 ranges, take into consideration the 
possibility of a complex phase where it may bring sizable contributions, obtained no significant dependency on the imaginary 
part of the new CKM factor. For the above mentioned quantities with a new family, deviations from the SM are promising, 
can be used as a probe of new physics. 
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I. INTRODUCTION 



Even if Standard Model (SM) is a successful theory, one should also check probable effects that may come from 
potential new physics. In the SM, since we do not have a clear theoretical argument to restrict number of generations 
to three, possibility of a new generation should not be ruled out until there is a certain evidence which order us to 
do so. This is especially true for rare B decays, which are very sensitive to generic expansions of the SM, due to 
their loop structure. We know from neutrino experiments that, for the mass of the extra generations there is a lower 
bound for the new generations (m„ 4 > 45 GeV) pj . Probable effects of extra generations was studied in many works 
The existing electroweak data on the Z-boson parameters, the W boson and the top quark masses excluded 
the existence of the new generations with all fermions heavier than the Z boson mass 16], nevertheless, the same 
data allows a few extra generations, if one allows neutral lcptons to have masses close to 50 GeV. In addition to this, 
recently observed neutrino oscillations requires an enlarged neutrino sector 0] . 

Generalizations of the SM can be used to introduce a new family, which was performed previously |l8j | . Using similar 
techniques, one can search fourth generation effects in B meson decays. The contributions from fourth generation 
to rare decays have been extensively studied 0, |2(], 152, H3| , where the measured decay rate has been used to 
put stringent constraints on the additional CKM matrix elements. In addition to B — > X s j , B — > X s l + £~ can be 
mentioned as one of the most promising areas in search of the fourth generation, via its indirect loop effects, to 
constrain V t * b Vf s The restrictions of the parameter space of nonstandard models based on LO analysis are 

not as sensitive as in the case of NLO analysis, hence a NLO analysis considering the possibility of a complex phase 
is important, which we plan to revise pfij . 

On the experimental side, the inclusive B — > XJ Jr £~ (with \fcp- > 0.2 GeV) decay with electron and muon modes 
combined {£ = e, /i) have been observed (Belle |27|). (BaBar J28|), 

B(B^X s £+£-) = (6.1 ± 1.4±^t) - (1) 
B(B -» X s £+£-) = (6.3 ± 1.6±i:|) • 10~ 6 . (2) 

They are in agreement with the SM B{B -> X s £+£-) SM = 4.2 ± 0.7 ■ 10~ 6 for the same cuts [29|. 

On the theoretical side, situation within and beyond the SM is well settled. A collective theoretical effort has led to 
the practical determination of B — > X s £ + £~ at the NNLO, which was completed recently, as a joint effort of different 
groups ( [30L I3lll32]| ) . and references therein. It is necessary to have precise calculations also in the extensions of the 
SM, which was performed for certain models. With the appearance of more accurate data we might be able to provide 
stringent constrains on the free parameters of the models beyond SM. From this respect, a NNLO analysis of the new 
generation is important. We study the contribution of the fourth generation in the rare B — ► X s £ + £~ decay at NNLO, 
to obtain experimentally measurable quantities which is expected to appear in the forthcoming years. 

The paper is organized as follows. In section 2, we present the necessary theoretical expressions for the 
B — > X s £ + f~decay in the SM with four generations. Section 3 is devoted to our conclusion. 

II. B -> X s e + r DECAY AND FOURTH GENERATION 

We use the framework of an effective low-energy theory, obtained by integrating out heavy degrees of freedoms, 
which in our case W-boson and top quark and an additional t' quark. Mass of the t' is at the order of fxw- I n this 
approximation the effective Hamiltonian relevant for the B — > X s £ + £~ decay reads [33| 

AC 10 

W eff =-- 7 IvZVt^GtoOiM , (3) 

* i—l 

where Gf is the Fermi coupling constant V is the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix, the 
the full set of th e op erators Oi (/i) and the corresponding expressions for the Wilson coefficients (7, (/x) in the SM can 
be found in Ref . |30| . 

In the model under consideration, the fourth generation is introduced in a similar way the three generations are 
introduced in the SM, no new operators appear and clearly the full operator set is exactly the same as in SM, which 
is a rough approximation. The fourth generation changes values of the Wilson coefficients Ci(n), i — 7, 8,9 and 10, 
via virtual exchange of the fourth generation up quark t' . With the definitions Xj = V* s Vjb, where j — u, c, t, t , the 
new physics Wilson coefficients can be written in the following form 

Cf G {pw) = -^-C l {llw)m t ~+ m t , , (4) 

M 
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where the last terms in these expression describes the contributions of the t' quark to the Wilson coefficients with the 
replacement of rrit with mt'- Notice that we use the definition \ t i — V t 1 s V t i b which is the element of the 4x4 Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, from now on '4G' will stand for sequential fourth generation model. In this model 
properties of the new t' quark are the same as ordinary t, except its mass and corresponding CKM couplings. A few 
comments are in order here: to obtain quantitative results we need the value of the fourth generation CKM matrix 
element V t * s Vt b which can be extracted i.e. from B — > X s j decay as a function of mass of the new top quark m' t . For 
this aim following 0113, we can use the fourth generation CKM factor Af in the range — 10~ 2 < At/ < 10~ 2 . In the 
numerical analysis, as a first step, Af is assumed real and expressions are obtained as a function of mass of the extra 
generation top quark m t >. It is interesting to notice that, if we assume Af can have imaginary parts, experimental 
values can also be satisfied |23T,l26| . Nevertheless, if we impose the unitarity condition of the CKM matrix we have 

V* s V ub + V* s V cb + V*V tb + V t 1 s V t , b = . (5) 

With the values of the CKM matrix elements in the SM ,34!], the sum of the first three terms in Eq. (5) is about 
7.6 x 1CP 2 , where the error in sum of first three terms is about ±0.6 x 1CP 2 . We assume the value of A t / is within 
this error range. 

What should not be ignored in constraining \ t > is that, when adding a fourth family the present constrains on 
the elements of CKM may get relaxed 34] . In order to have a clear picture of X t >, CKM matrix elements should 
be calculated with the possibility of a new family, using present experiments that constitutes the CKM. From this 
respect we do not have to exclude certain regions that violate unitarity of the present CKM, but take it in the ranges 
-lO" 2 < A t < < 10- 2 and -10" 3 < X t , < ICT 3 . 



A. Differential Decay Width 



Since extended models are very sensitive to NNLO corrections, we used the NNLO expression for the branching 
ratio of the radiative decay B — > X s £ + £~, which has been presented in Refs. |2^,|3^]. In the NNLO approximation, 
the invariant dilepton mass distribution for the inclusive decay B — > X s £ + £~ can be written as 



dT(b -> X s £ + n (ae m \ 2 G' F ml pole \V* s V tb \ 

ds V 4tt J 48tt 3 [ S) 



x (1 + 2s) 



a 



eff 



^cff 



+ 4 (1 + 2/J) 



off 



cff^cffs 



12Re C? L C\ 



(6) 



where s = mj + ^/ml pole with (£ = e or fi). In the SM the effective Wilson coefficients Cq and Cfg are given 

by [3(1133 and can be obtained from Eqs.JBEl and [TDJ, by setting 4G — > 0. Following the lines of A.Ali (2^ with 
the assumption that only the lowest non-trivial order of these Wilson coefficients get modified by new physics, which 

means that C^}\nw), C^d^w), C^d^w) and C^{^w) get modified, the shifts of the Wilson coefficients at fiw 
can be written as 



CiM — > C t ( P w) + j 1 Cfiiiw) ■ 
These shift at the matching scale are resulted in the modifications of the effective Wilson coefficients, 



(7) 



Cf 



1 + ^iM W7 (j) ) {Al + A 77 C$ G ( m ) + A 7S Ct G {n w )) 
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c: 



°io — 



1 + ^1 W9 (|)J (a, + T 9 h(m 2 c , S) + U 9 h(l, s) + W 9 h(0, s) + C 9 4G (^)) 



7T 



The numerical values for the parameters A 77 , A 7 %, A$g , which incorporate the effects from the running, can be found 
in the same reference |29). for the functions /i(m 2 , s) and tog(s), they are given in Ref. [13, while lo 7 (s) and F^(s) 
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FIG. 1: Branching ratio B 



B->x s e+t~ 



[10~ 6 ] as a function of s e [0.05, 0.25](see Eq.JTTJ)- The four thick lines show the NNLL 



prediction for rati = 200, 300, 400 and 500 with increasing thickness respectively and the SM prediction is the thin line. The 



figures are obtained at the scale /i = 5.0 GeV. For the figure at the Left: At 



-10 , Right: At, = 10~ 
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FIG. 2: The same as Fig^with the choices, For the figure at the Left: A t , = —10 3 , Right: A t , = 10 



can be seen in Ref. 1.331. In order to remove the large uncertainty coming from m& terms it is customary to use the 
following expression 29] 



B 



(s) 



®t7 Xc ™ dT(B -> X s £+£-) 
T(B X c ev) d~s 



(11) 



which can be called as branching ratio. The explicit expression for the semi-leptonic decay width can be found in Ref. 
[3fj|. The branching ratio with 4G is presented in Figs. for the choice of the scale /i = 5 GeV. 

In the figures related with dilepton invariant mass distribution we used the low region s € [0.05, 0.25] where peaks 
stemming from cc resonances are expected to be small. During the calculations we take B^ > T* Xcel ' — 0.1045. 



B. Forward-Backward asymmetry 



The double differential decay width d 2 T(b — > X s £ + £ )/(dsdz), (z — cos(0)) is expressed as 



iymn 
31] 



We investigate both, the so-called normalized and the unnormalizcd forward-backward asymmetry with 4G model. 

K s £+£-)/(dsdz), (; 

-i2 ^5 n/*T/. i 2 



d 2 T(b -» X s £+£-) (a em \2 G 2 F m b bpole \V* s V tb 



ds dz 



( U em y '"-b,polc I v ts v tb\ _ 2 

UttJ 48tt3 ( S) 



4 



:[(l-z 2 ) + .s(l + z 2 )] 



a 



off 



+;[(! + z 2 ) + s(l-z 2 )} Cf (l+^£/ 77 ( s » 



c 



2a s 



f9s(s,Z) 
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-3 s z Re(C 9 eff C 1 f*) f 1 + H2i / 910 (|) 

7T 
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7T 



2a s 
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2a s 



(12) 



where 9 is the angle between the momenta of the b quark and the £ + , measured in the rest frame of the lepton pair. 
The functions fgg(s,z), fn(s,z), /gm(s), f-jg{s,z) and /Vio(s) are the analogues of Wgg(s), (jJtt{s) and W7g(s) which 
can be found in the same reference l31l 




FIG. 3: Unnormalized forward-backward asymmetry Ayb. [10 ] as a function of s £ [0,0.25] (see Eq.JTJl). The four thick 
lines show the NNLL prediction for rati = 200, 300, 400 and 500 with increasing thickness respectively and the SM prediction 
is the thin line. The figures are obtained at the scale \i — 5.0 GeV. For the figure at the Left: At, = — 10~ 2 , Right: At, = 10 -2 . 
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FIG. 4: The same as FigEJwith the choices: Left: A t , = -10~ 3 , Right: At, = 10" 3 



The unnormalized version of forward-backward asymmetry, ^4fb(s) is defined as 



Afb(s) 



J-r 



1 d 2 r(b^X s 



ds dz 



■ sgn(z) dz ^ 



T{B -» X c eD e ) 



B 



(13) 
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FIG. 5: Normalized forward-backward asymmetry Afb [10 -6 ] as a function of s £ [0,0.25] (see Eg. 1141 '). The four thick lines 
show the NNLL prediction for mti = 200, 300, 400 and 500 with increasing thickness respectively and the SM prediction is the 
thin line. The figures are obtained at the scale = 5.0 GeV . For the figure at the Left: At/ = — 10~ 2 , Right: At/ = 10 -2 . 




while the definition of the normalized forward-backward asymmetry Ayb(s) reads 



1 <*\ /-i dsdz ' sgn(z)dz 

Afb ^ s > = — ri dT^x.t+rr ^: — ■ [ > 

J — X dsdz 

The position of the zero of the Afb{so) = is very sensitive to 4G effects as it is seen in the figures 1)315(1 . However 
as 4G parameter Ac decreases expectations of the new model are getting closer to SM values which can be inferred 
from Figs. 1(416)1 



C. Integrated Branching Ratio 



By suitable choice of integration limits over s one can obtain integrated branching ratio in accordance with the 
experiment for e and u, which is already performed, hence we use the integrated branching ratio expression which 
has the following form 29] : 



B[B -> X s £+r) = I(T 6 x L + a 2 \A 7 ot \ 2 + a 3 (|C 9 4G | 2 



\<\ 2 ) 



+a 4 Re^ * ReG 
+a 7 ImA^ - 



•4G 



■0 8 1m AZ* 



ReCf 3 



a 9 ImCg G 



- a 6 Re A\' 
oio Re C 10 



tot 



(15) 
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FIG. 7: Integrated Branching ratio B(B — > X s £ + £~) [10~ 6 ] as a function of \ t , for £ — e (see Eq. 11511 '). In the left figure 
At, G [— 10 -2 , 10~ 2 ]. For the figure at the right A t , 6 [— 1CP 3 , 10~ 3 ]. In the figures straight lines shows the SM allowed region. 
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FIG. 8: Integrated Branching ratio B(B -> X s £ + £~) [10~ 6 ] as a function of At, for I = fi. In the left figure A t , e [-10 -2 , 1(T 2 ]. 
For the figure at the right A t , £ [— 10~ 3 , 10 -3 ]. In the figures straight lines show the SM region. 

where the numerical value of the coefficients ai are given in Table Q]for £ = e, /i. For the integrated branching 
ratios we refer to Figs. (|7l8|l of electron and muon respectively. 
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ai 


a 2 




C14 


as 


a@ 


CL7 


as 


ag 


aio 


e 


1.9927 


6.9357 


0.0640 


0.5285 


0.6574 


0.2673 


-0.0586 


0.4884 


0.0095 


-0.5288 


/' 


2.3779 


6.9295 


0.0753 


0.6005 


0.7461 


0.5955 


-0.0600 


0.5828 


0.0102 


-0.6225 



TABLE I: Numerical values of the coefficients ai (evaluated at )j> b = 5 GeV) for the decays B — > X s £ + £ (£ = e, p), taken from 
Ref. HJ. 



III. DISCUSSION 

In the sequential fourth generation model, there are basically two free parameters, mass of new generations and 
CKM factors which can have imaginary phases. As a worst scenario, we decompose X t/ = Re[Xt/] + 1 x Im[\ tl ] and 
choose the range ^rxffr — 10~ 2 ; we checked the effect of this choice and observe that contribution from the imaginary 
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part can be neglected for all of the kinematical observables. Naturally, these quantities should be fixed by respecting 
experiment. Besides, constraints for CKM values should be updated by noting that existance of a new generation can 
relax the matrix elements of CKM3 X 3, when it is accepted as a sub-matrix of CKM/ty4. 

Since scale dependency of NNLO calculations of B — > X s £ + £~ are not very high |3l|, during the calculations we 
set the scale /x = 5 GeV, use the main input parameters as follows, 

a em = 1/133 , a s (m z ) = 0.119 , Gf = 1.16639 x 10~ 5 GeV~ 2 , m w — 80.33 GeV , 
rrib = 4.8GeV,TO t = 176 GeV, m c = 1.4 GeV, Wolfenstcin parameters: 
A = 0.75, A = 0.221, p= 0.4, r\ = 0.2. (16) 

Effects of new physics on kinemaical observables can be summerized as follows: 

• Differential decay width Q B ~* X ^ + ^ [ s presented in figures Fig. (|l!2f) . where it is shown that SM prediction can 
be strongly enhanced with a new quark for the choice X tf < 0. It is also possible to supress the decay width for 
positive solutions of Xt> which is not favored. 

• Forward-Backward asymmetry is also very sensitive to 4G effects, especially for the choice At/ = 10~ 2 . As it is 
seen in Figs.JSEJ), as the mass of mc increases it is even possible to have positive values for Ap s(0) which is in 
contradiction with SM, but natural in extended models. Once the experimental results related with this quantity 
is obtained, it will be a keen test of fourth generation model. Deviations from the point s=0 are detectable as it 
is seen in Fig.Q for the choice of Xn G [— 10~ 3 , 10 -3 ], whereas for the same region we see almost no dependence 
on the normalized forward-backward asymmetry in Fig.@. While Standard Model states the central value 
^fb LO (°) = -(2-30 ± 0.10) x 10~ 6 , 4G predictions cover the range ^fb' NNLO (0) G [-6, 1] x 10~ 6 for the choices 
Xt> = — 10 -2 , 10 -2 respectively. For the point where forward-backward asymmetry vanishes Standard Model 
result is s^ NLO = 0.162 ± 0.002 however 4G predictions are roughly ^g.nnlo g [0.13,0.18]. 

• Integrated branching ratios Figs. (|7I8|) strongly depends on the new physics parameters Xv and mp, therefore it 
is possible to restrict them by respecting experiments. As it can be deduced from the figures when 4G effects 
are switched off our calculations are lying on the SM ground within error bars [29| . Similar to branching ratio 
for integrated branching ratios enhancement comes from negative choices of Af which favors smaller values for 
a sm,nnlo (0) = _(2. 30 ±0.10) x 10- 6 . 

To summarize, in this work we present the predictions of the sequential fourth generation model for experimentally 
measurable quantities related with B — > X s £ + £~ decay which is expected to emerge in the near future thanks to 
running B factories. These predictions differ from SM in certain regions, hence can be used, to differentiate the 
existence of the fourth family or to put stringent constrains on the free parameters of the model, if it exists. 
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